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The directed flow of charged hadron and identified particles has been studied in the framework 
of a multi-phase transport (AMPT) model, for ^*"'Au+^^''Au collisions at =200, 130, 62.4, 

39, 17.2 and 9.2 GeV. The rapidity, centrality and energy dependence of directed flow for charged 
particles over a wide rapidity range are presented. AMPT model gives the correct slope, as 

well as its trend as a function of energy, while it underestimates the magnitude. Within the AMPT 
model, the proton v\ slope is found to change its sign when the energy increases to 130 GeV - a 
feature that is consistent with "anti-flow". Hadronic re-scattering is found having little effect on 
Hi at top RHIC energies. These studies can help us to understand the collective dynamics at early 
times in relativistic heavy-ion collisions, and they can also be served as references for the RHIC 
Beam Energy Scan program. 

PACS numbers: 25.75.Ld, 25.75.Nq, 25.75.Dw 



I. INTRODUCTION 



Anisotropic flow is one of the key observables in charac- 
terizing properties of the dense and hot medium created 
in the relativistic heavy-ion collisions [T|. It is quanti- 
fied by Fourier coefficients when expanding particle's az- 
imuthal distribution with respect to the reaction planeQ: 
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(Pp 2-K pxdpTdy 
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where (/> denotes the angle between the particle's az- 
imuthal angle in momentum space and the reaction plane 
angle. The sine terms in Fourier expansions vanish due 
to the reflection symmetry with respect to the reaction 
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plane. The various coefficients in this expansion can be 
defined as: 



(cos 71(7 



(2) 



The first and the second coefficients are named as di- 
rected flow (vi) and elliptic flow(i;2), respectively, and 
they play important roles in describing the collective ex- 
pansion in azimuthal space. Elliptic flow is produced by 
the conversion of the initial coordinate-space anisotropy 
into momentum-space anisotropy, due to the developed 
large in-plane pressure gradient. Elliptic flow depends 
strongly on the re-scattering of the system constituents, 
thus it is sensitive to the degree of thermalizatioufs*! of 
the system at early time. Directed fiow, which is the fo- 
cus of this study, describes the "side splash" of particles 
away from mid-rapidity [S], and it probes the dynamics of 
the system in the longitudinal direction. Since directed 
flow is generated very early, it brings information from 
the foremost early collective motion of the system. The 
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shapes of directed flow, in particular those for identified 
particles, are of special interest because they are sensi- 
tive to the equation of state (EOS) and may carry a phase 
transition sie nalf^. 

The study of energy dependence of directed flow has 
implications in many aspects. Firstly, because directed 
flow has a unique, pre-cquilibrium origin, it is expected 
to behave differently than other soft observables which 
show an "entropy-driven" multiplicity scaling Q. It has 
been shown by STARf?'] that at top RHIC energies, di- 
rected flow is independent of system size, while it has 
an energy dependence. For a comprehensive study of 
the subject, it is necessary to extend the study of the 
energy dependence of directed flow in a wider energy 
range. Secondly, experiments at RHIC (PHENIX and 
STAR) have planned to look for the existence of the QCD 
phase boundary and the possible critical point by collid- 
ing heavy ions at various incident beam energies [a. [ol. [lOj . 
A non-monotonic dependence of variables on ^/snn and 
an increase in event-by-event fluctuations should become 
apparent near the critical point [sj . Directed flow is gener- 
ated during the nuclear passage time {2R/j ~ 0.1/m/c) 
and it probes the onset of bulk collective dynamics in the 
earlier stage of the collision. As a suggested signature of 
a first order phase transition Q, directed fiow is sensitive 
to the creation of the critical point and it plays an im- 
portant role in the proposed beam energy scan program. 

In this paper, directed flow from the AMPT model for 
6 energies are presented. They are, 9.2, 17.3, 39, 62.4, 
130 and 200 GeV. The comparison with the measure- 
ments from STAR and PHOBOS are made at top ener- 
gies. The particle type dependence over a wide rapidity 
range is discussed. This study will deepen our under- 
standing about the energy dependence of directed flow, 
and it can be also served as a valuable reference for the 
RHIC Beam Energy Scan program. 
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FIG. 1: Rapidity dependence of vi for charged particles 
in the AMPT model compared with STAR and PHOBOS 
data (plotted as a function of rj) in the Au-|-Au collisions at 
ysiTzv = 200GeV. The dashed lines showed AMPT result 
from different centrality: 0-30% (black), 30%-60%(red), 60%- 
80%(blue). 



model[l^ with modifications and extensions. As sug- 
gested in Ref . [l9| . the parton cross section is chosen as 
3 mb in our analysis. All the errors presented here are 
statistical only. 



II. THE AMPT MODEL 

The AMPT model consists of four main 
components[ll| : initial conditions, partonic interac- 
tions, conversion from partonic matter to hadronic 
matter, and hadronic interactions. The initial con- 
ditions, which include the spatial and momentum 
distributions of the mini-jet partons and soft string 
excitations, are obtained from the HIJING model p^. 
The scatterings among partons are modeled by Zhang's 
parton cascade(ZPC)jl4j, which includes two-body scat- 
terings with cross sections from pQCD with screening 
masses. In the default AMPT model[lB|, partons are 
recombined with their parent strings when they stop 
interacting, and the resulting strings fragment into 
hadrons according to the Lund string fragmentation 
model [ill . In the AMPT model with string melting [13, 
quark coalescence is used instead to combine partons 
into hadrons. The dynamics of the subsequent hadronic 
matter is described by the ART (a relativistic transport) 



III. ANALYSIS AND RESULTS 

In Fig. [U the directed flow of charged particles from 
AMPT is shown as a function of rapidity, for collision en- 
ergies of 200, 130, 62.4, 39, 17.2 and 9.2GeV. The central- 
ity is divided into three bins, namely, 0-30%, 30%-6G% 
and 60%-80%, based on the impact parameter (6) dis- 
tribution. The calculations with string melting scenario 
is used for high energies (200, 130, 62.4, 39GeV) while 
for low energies (17.2 and 9.2GeV) calculations are per- 
formed with default scenario. The reason for such choice 
is because that, it is argued [l^, [H, [llj that the string 
melting should be used to explain flow around midrapid- 
ity at top RHIC energies, and default setting describes 
data at 9.2 GeV the best. The energy density in the col- 
lisions at the RHIC top energies is mush higher than the 
critical density for the QCD phase transition. More dis- 
cussion on different AMPT configurations can be found 
later in this paper. All results are obtained by integrating 
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FIG. 3: Cliarged particles vi{y) from centrality 10%-70% 
for 9.2GeV(upper left panel), 17. 3GeV(upper right panel) and 
39GeV(down left panel). The dashed lines show three AMPT 
versions : string melting scenario(black), default scenario 
with high-NTMAX(red) and low-NTMAX(blue). Experimen- 
tal data points from STAR are plotted as a function of rj. 



FIG. 2: Proton(solid lines) and Pion(dashed lines) vi{y) from 
AMPT at centrality 10%-70%. 



over transverse momentum (pr) up to 4.0 GeV/c. Exper- 
imental results from STARp, ^ and PHOBOSlH] are 
also shown for comparison. The charged hadron vi mea- 
sured by the PHOBOS experiment is for 0-40% central 
collisions, and the results measured by STAR experiment 
are for centrality 30%-60% at 200GeV, and centrality 0- 
60% at 9.2GeV. In general, AMPT gives larger vi at 
low energies than at high energies, the same trend has 
been seen in data. At top RHIC energies, AMPT un- 
derestimated vi , due to the turn-off of mean-field poten- 
tials in ART when implemented in AMPT to describe the 
hadronic scattering However, in the rapidity range 
of [-2.0,2.0], the shape of vi between AMPT calculations 
and experimental data are in good agreement — this can 
be seen by scaling experimental results with a factor of 
0.25. 

The particle type dependence of directed flow is shown 
in Fig. O The different sign of vi between pions and 
protons at low energies can be understood as nucleon 
shadowing and baryon stopping[23. l25l |. In general the 
magnitude of the vi slope at midrapidity decreases with 
increasing energy. This effect is most profound for pro- 
tons, for which the slope keeps decreasing and when the 
energy is high enough, it changes its sign and protons be- 
gin to flow together with pions. This is consistent with 
the "anti-flow" scenario [26|, in which the "bounce-off" 
motion and transverse expansion of nucleons compete 
with each other around midrapidity, and when the trans- 
verse expansion is strong enough (e.g., at top RHIC en- 
ergies) , it overcomes the "bounce-off" motion and causes 



protons to change their sign of directed flow. 

To illustrate the effect on vi due to different configu- 
rations in AMPT, in Fig. [3] we present the directed flow 
of charged hadrons in low energy collisions, from AMPT 
calculations with the string melting scenario and the de- 
fault scenario. The similar study for higher energies can 
be found in (loj . The calculation with string melting 
yields the smallest vi slope around mid-rapidity and is 
close to data. Two different default scenarios are also 
studied: one is calculated with NTMAX=2500 (high- 
NTMAX), and the other, NTMAX=150 (low-NTMAX). 
NTMAX stands for the number of time-steps for the 
hadron cascade (see detail in paper [ll|). A large NT- 
MAX means a thoroughly developed hadron cascade, as 
0.2fm/c*NTMAX is the termination time, in the center 
of mass frame, of the hadron cascade in AMPT model. 
The comparison, for low energies, of vi calculated be- 
tween low-NTMAX and high-NTMAX indicates that vi 
can change its sign at large rapidity if the time for the 
hardonic cascade is long enough. In default AMPT, the 
NTMAX has to be much larger than 150 in order to de- 
scribe vi at large rapidity. The disagreement between the 
experimental data and the calculation made with high- 
NTMAX is mostly due to the lack of the mean-field in the 
hadron cascade in AMPT, which is a considerable effect 
at low energies when the nuclei passage time is not neg- 
ligible (compared to that at high energies). The AMPT 
calculation with high-NTMAX at high energy has been 
presented in [l2j. In this paper, we address the compar- 
ison around midrapidity only, and results presented in 
this paper are made with low-NTMAX unless otherwise 
specified. 

The energy dependence of charged particle directed 
fiow, calculated with the AMPT model, is shown in 
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FIG. 4: Charged hadrons' slope dvi/dy' in the mid-rapidity 
\y'\ < 0.5 as a function of incident-energy. The data are taken 
from STAR(stars) and PHOBOS(squares) and scaled by a 
factor 0.25. The AMPT calculatons with string melting before 
ART are depicted with open circles and after hadron cascade 
are depicted with full circles. The open triangles depict the 
default AMPT calculaitons before ART and the full triangles 
depict after hadron cascade. 



Fig. m Experimental data are also shown for compar- 
ison. The centrality for which the calculation is per- 
formed is 10%-70%. The centrality for PHOBOS data 
from different energies is 0-40% while the centrality se- 
lection for STAR data are 0-60% for 9.2 GeV, 10%-70% 
for 62.4 GeV, and 30%-60% for 200 GeV. To obtain the 
integrated vi , one needs to fold in the spectra at different 
energies, which brings in an additional layer of systemat- 
ics. Thus instead, we present the slope of vi{y) around 
mid-rapidity {\y'\ < 0.5) extracted from the normalized 
in' = y/Ubeam) rapidity distribution, where ybeam is the 
beam rapidity. For the energy range that string melting 
is used (39 GeV and above), all the AMPT calculations 
underestimate the experimental data, however, they pre- 
dict the right trend of the energy dependence. For the 
low energies at 9.2 GeV, calculations with string melting 
did a poor job, the calculation with the default AMPT 



improves the result in the right direction yet is still not 
be able to explain the data. The hadron re-scattering 
effect on directed flow vi can be seen by switching off the 
hadron cascade in the AMPT calculation. Comparing the 
difference between the result with hadron cascade (open 
symbols) and without (solid symbols), it is found that the 
hadronic cascade has a significant effect for low energy 
results but little for that of high energies. This can be 
understood as that, when the energy is high enough, the 
hadron re-scattering become less important due to the 
presence of strong collective motion built up beforehand. 

IV. SUMMARY AND CONCLUSIONS 

In this paper, vi values calculated from the AMPT 
model for different energies are discussed. It is found 
that the AMPT model gives the right shape of vi versus 
y while underestimating the magnitude, possibly due to 
the lack of mean-field in its hadron cascade. In AMPT, 
the proton vi slope changes its sign when the energy in- 
creases to 130 GeV and begins to have the same sign 
as that of pions, as expected in the "anti-flow" scenario. 
The effect on vi due to string melting, low-NTMAX and 
high-NTMAX are illustrated. The energy dependence of 
the vi slope at midrapidity is compared to experimental 
data, and AMPT can describe the trend of energy de- 
pendence while missing the magnitude by a fraction of 
75%. Hadronic rescattering is found to be less important 
at high energies as the strong collective motion becomes 
to be the dominant dynamics. 
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I. INTRODUCTION 



Anisotropic flow is one of the key observables in charac- 
terizing properties of the dense and hot medium created 
in the relativistic heavy-ion collisions [T|. It is quanti- 
fied by Fourier coefficients when expanding particle's az- 
imuthal distribution with respect to the reaction planeQ: 
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angle. The sine terms in Fourier expansions vanish due 
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plane. The various coefficients in this expansion can be 
defined as: 
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The first and the second coefficients are named as di- 
rected flow (vi) and elliptic flow(i;2), respectively, and 
they play important roles in describing the collective ex- 
pansion in azimuthal space. Elliptic flow is produced by 
the conversion of the initial coordinate-space anisotropy 
into momentum-space anisotropy, due to the developed 
large in-plane pressure gradient. Elliptic flow depends 
strongly on the re-scattering of the system constituents, 
thus it is sensitive to the degree of thermalizatioufs*! of 
the system at early time. Directed fiow, which is the fo- 
cus of this study, describes the "side splash" of particles 
away from mid-rapidity [S], and it probes the dynamics of 
the system in the longitudinal direction. Since directed 
flow is generated very early, it brings information from 
the foremost early collective motion of the system. The 
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shapes of directed flow, in particular those for identified 
particles, are of special interest because they are sensi- 
tive to the equation of state (EOS) and may carry a phase 
transition sie nalf^. 

The study of energy dependence of directed flow has 
implications in many aspects. Firstly, because directed 
flow has a unique, pre-cquilibrium origin, it is expected 
to behave differently than other soft observables which 
show an "entropy-driven" multiplicity scaling Q. It has 
been shown by STARf?'] that at top RHIC energies, di- 
rected flow is independent of system size, while it has 
an energy dependence. For a comprehensive study of 
the subject, it is necessary to extend the study of the 
energy dependence of directed flow in a wider energy 
range. Secondly, experiments at RHIC (PHENIX and 
STAR) have planned to look for the existence of the QCD 
phase boundary and the possible critical point by collid- 
ing heavy ions at various incident beam energies [a. [ol. [lOj . 
A non-monotonic dependence of variables on ^/snn and 
an increase in event-by-event fluctuations should become 
apparent near the critical point [sj . Directed flow is gener- 
ated during the nuclear passage time {2R/j ~ 0.1/m/c) 
and it probes the onset of bulk collective dynamics in the 
earlier stage of the collision. As a suggested signature of 
a first order phase transition Q, directed fiow is sensitive 
to the creation of the critical point and it plays an im- 
portant role in the proposed beam energy scan program. 

In this paper, directed flow from the AMPT model for 
6 energies are presented. They are, 9.2, 17.3, 39, 62.4, 
130 and 200 GeV. The comparison with the measure- 
ments from STAR and PHOBOS are made at top ener- 
gies. The particle type dependence over a wide rapidity 
range is discussed. This study will deepen our under- 
standing about the energy dependence of directed flow, 
and it can be also served as a valuable reference for the 
RHIC Beam Energy Scan program. 
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FIG. 1: Rapidity dependence of vi for charged particles 
in the AMPT model compared with STAR and PHOBOS 
data (plotted as a function of rj) in the Au-|-Au collisions at 
ysiTzv = 200GeV. The dashed lines showed AMPT result 
from different centrality: 0-30% (black), 30%-60%(red), 60%- 
80%(blue). 



model[l^ with modifications and extensions. As sug- 
gested in Ref . [l9| . the parton cross section is chosen as 
3 mb in our analysis. All the errors presented here are 
statistical only. 



II. THE AMPT MODEL 

The AMPT model consists of four main 
components[ll| : initial conditions, partonic interac- 
tions, conversion from partonic matter to hadronic 
matter, and hadronic interactions. The initial con- 
ditions, which include the spatial and momentum 
distributions of the mini-jet partons and soft string 
excitations, are obtained from the HIJING model p^. 
The scatterings among partons are modeled by Zhang's 
parton cascade(ZPC)jl4j, which includes two-body scat- 
terings with cross sections from pQCD with screening 
masses. In the default AMPT model[lB|, partons are 
recombined with their parent strings when they stop 
interacting, and the resulting strings fragment into 
hadrons according to the Lund string fragmentation 
model [ill . In the AMPT model with string melting [13, 
quark coalescence is used instead to combine partons 
into hadrons. The dynamics of the subsequent hadronic 
matter is described by the ART (a relativistic transport) 



III. ANALYSIS AND RESULTS 

In Fig. [U the directed flow of charged particles from 
AMPT is shown as a function of rapidity, for collision en- 
ergies of 200, 130, 62.4, 39, 17.2 and 9.2GeV. The central- 
ity is divided into three bins, namely, 0-30%, 30%-6G% 
and 60%-80%, based on the impact parameter (6) dis- 
tribution. The calculations with string melting scenario 
is used for high energies (200, 130, 62.4, 39GeV) while 
for low energies (17.2 and 9.2GeV) calculations are per- 
formed with default scenario. The reason for such choice 
is because that, it is argued [l^, [H, [llj that the string 
melting should be used to explain flow around midrapid- 
ity at top RHIC energies, and default setting describes 
data at 9.2 GeV the best. The energy density in the col- 
lisions at the RHIC top energies is mush higher than the 
critical density for the QCD phase transition. More dis- 
cussion on different AMPT configurations can be found 
later in this paper. All results are obtained by integrating 



3 



AMPT Melting 
— 71 + 71" 
proton 



-5 
5 



200 GeV:: 



62.4 GeV:: 



20 

-20 



AMPf belault' 
proton 



17.3GeV: 



10%-70% p^<4.0 



130 GeV: 



39 GeV 



9.2 GeV 



10 

5 

-5 

g-10 
> 5 

-5 







9.2GeV 



\ 

\ 



39 GeV 



17.3GeV 



AMPT10%-70% p.^<4.0 
melting 

default high-NTMAX 

default low-NTMAX 

□ STAR X 0.25 



-4 -2 2 4 , .-4 -2 2 4 

y (11) 



FIG. 3: Cliarged particles vi{y) from centrality 10%-70% 
for 9.2GeV(upper left panel), 17. 3GeV(upper right panel) and 
39GeV(down left panel). The dashed lines show three AMPT 
versions : string melting scenario(black), default scenario 
with high-NTMAX(red) and low-NTMAX(blue). Experimen- 
tal data points from STAR are plotted as a function of rj. 



FIG. 2: Proton(solid lines) and Pion(dashed lines) vi{y) from 
AMPT at centrality 10%-70%. 



over transverse momentum (pr) up to 4.0 GeV/c. Exper- 
imental results from STARp, ^ and PHOBOSlH] are 
also shown for comparison. The charged hadron vi mea- 
sured by the PHOBOS experiment is for 0-40% central 
collisions, and the results measured by STAR experiment 
are for centrality 30%-60% at 200GeV, and centrality 0- 
60% at 9.2GeV. In general, AMPT gives larger vi at 
low energies than at high energies, the same trend has 
been seen in data. At top RHIC energies, AMPT un- 
derestimated vi , due to the turn-off of mean-field poten- 
tials in ART when implemented in AMPT to describe the 
hadronic scattering However, in the rapidity range 
of [-2.0,2.0], the shape of vi between AMPT calculations 
and experimental data are in good agreement — this can 
be seen by scaling experimental results with a factor of 
0.25. 

The particle type dependence of directed flow is shown 
in Fig. O The different sign of vi between pions and 
protons at low energies can be understood as nucleon 
shadowing and baryon stopping[23. l25l |. In general the 
magnitude of the vi slope at midrapidity decreases with 
increasing energy. This effect is most profound for pro- 
tons, for which the slope keeps decreasing and when the 
energy is high enough, it changes its sign and protons be- 
gin to flow together with pions. This is consistent with 
the "anti-flow" scenario [26|, in which the "bounce-off" 
motion and transverse expansion of nucleons compete 
with each other around midrapidity, and when the trans- 
verse expansion is strong enough (e.g., at top RHIC en- 
ergies) , it overcomes the "bounce-off" motion and causes 



protons to change their sign of directed flow. 

To illustrate the effect on vi due to different configu- 
rations in AMPT, in Fig. [3] we present the directed flow 
of charged hadrons in low energy collisions, from AMPT 
calculations with the string melting scenario and the de- 
fault scenario. The similar study for higher energies can 
be found in (loj . The calculation with string melting 
yields the smallest vi slope around mid-rapidity and is 
close to data. Two different default scenarios are also 
studied: one is calculated with NTMAX=2500 (high- 
NTMAX), and the other, NTMAX=150 (low-NTMAX). 
NTMAX stands for the number of time-steps for the 
hadron cascade (see detail in paper [ll|). A large NT- 
MAX means a thoroughly developed hadron cascade, as 
0.2fm/c*NTMAX is the termination time, in the center 
of mass frame, of the hadron cascade in AMPT model. 
The comparison, for low energies, of vi calculated be- 
tween low-NTMAX and high-NTMAX indicates that vi 
can change its sign at large rapidity if the time for the 
hardonic cascade is long enough. In default AMPT, the 
NTMAX has to be much larger than 150 in order to de- 
scribe vi at large rapidity. The disagreement between the 
experimental data and the calculation made with high- 
NTMAX is mostly due to the lack of the mean-field in the 
hadron cascade in AMPT, which is a considerable effect 
at low energies when the nuclei passage time is not neg- 
ligible (compared to that at high energies). The AMPT 
calculation with high-NTMAX at high energy has been 
presented in [l2j. In this paper, we address the compar- 
ison around midrapidity only, and results presented in 
this paper are made with low-NTMAX unless otherwise 
specified. 

The energy dependence of charged particle directed 
fiow, calculated with the AMPT model, is shown in 
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FIG. 4: Charged hadrons' slope dvi/dy' in the mid-rapidity 
\y'\ < 0.5 as a function of incident-energy. The data are taken 
from STAR(stars) and PHOBOS(squares) and scaled by a 
factor 0.25. The AMPT calculatons with string melting before 
ART are depicted with open circles and after hadron cascade 
are depicted with full circles. The open triangles depict the 
default AMPT calculaitons before ART and the full triangles 
depict after hadron cascade. 



Fig. m Experimental data are also shown for compar- 
ison. The centrality for which the calculation is per- 
formed is 10%-70%. The centrality for PHOBOS data 
from different energies is 0-40% while the centrality se- 
lection for STAR data are 0-60% for 9.2 GeV, 10%-70% 
for 62.4 GeV, and 30%-60% for 200 GeV. To obtain the 
integrated vi , one needs to fold in the spectra at different 
energies, which brings in an additional layer of systemat- 
ics. Thus instead, we present the slope of vi{y) around 
mid-rapidity {\y'\ < 0.5) extracted from the normalized 
in' = y/Ubeam) rapidity distribution, where ybeam is the 
beam rapidity. For the energy range that string melting 
is used (39 GeV and above), all the AMPT calculations 
underestimate the experimental data, however, they pre- 
dict the right trend of the energy dependence. For the 
low energies at 9.2 GeV, calculations with string melting 
did a poor job, the calculation with the default AMPT 



improves the result in the right direction yet is still not 
be able to explain the data. The hadron re-scattering 
effect on directed flow vi can be seen by switching off the 
hadron cascade in the AMPT calculation. Comparing the 
difference between the result with hadron cascade (open 
symbols) and without (solid symbols), it is found that the 
hadronic cascade has a significant effect for low energy 
results but little for that of high energies. This can be 
understood as that, when the energy is high enough, the 
hadron re-scattering become less important due to the 
presence of strong collective motion built up beforehand. 

IV. SUMMARY AND CONCLUSIONS 

In this paper, vi values calculated from the AMPT 
model for different energies are discussed. It is found 
that the AMPT model gives the right shape of vi versus 
y while underestimating the magnitude, possibly due to 
the lack of mean-field in its hadron cascade. In AMPT, 
the proton vi slope changes its sign when the energy in- 
creases to 130 GeV and begins to have the same sign 
as that of pions, as expected in the "anti-flow" scenario. 
The effect on vi due to string melting, low-NTMAX and 
high-NTMAX are illustrated. The energy dependence of 
the vi slope at midrapidity is compared to experimental 
data, and AMPT can describe the trend of energy de- 
pendence while missing the magnitude by a fraction of 
75%. Hadronic rescattering is found to be less important 
at high energies as the strong collective motion becomes 
to be the dominant dynamics. 
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